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1. Introduction 

Colour reconnection (also referred to as 'rearrangement' or 'recoupling') in W+W~ 
decays has been the subject of many studies (e.g. ||, ^) and at present there is 
agreement that observable effects of interference between the colour singlets in the 
perturbative phase are expected to be small. In contrast, significant interference 
in the hadronisation process appears a viable prospect but, with our current lack 
of knowledge of non-perturbative qcd, such interference can only be estimated 
in the context of specific models H, ^, ||, ||, |^, |^. In the studies described 
below, experimentally accessible features of these models^ are investigated, paying 
particular attention to the bias introduced to a typical measurement of Afyv by direct 
reconstruction of the decay products. 

Throughout this section reconnection effects were studied using: pythia 5.722, 
type I and type II superconductor models (with the string overlap integral in type 
I case characterised by p = 0.9) [|[ ||]; ariadne 4.08 allowing reconnection between 
the two W bosons; and herwig 5.9, in both its default reconnection model and also 
a 'colour octet' variant in which merging of partons to form clusters was performed 
on a nearest neighbour basis|. In all cases, the tuning of the models was as used in 
reference M. 



2. Inclusive charged multiplicity 

It has been suggested [Q, ^ that simple observable quantities such as the charged 
multiplicity in restricted rapidity intervals may be sensitive to the effects of colour 
reconnection. More recently [g| it was suggested that the effect on the inclusive 
charged multiplicity itself may be larger than previously considered and that the 
mean hadronic multiplicity in W^W^— >qqqq events, (N^^), may be as much as 10% 
smaller than twice the hadronic multiplicity in W+W~^qq£F£ events, {N^i^^'^)- It 
was also reported during this workshop that the effects of Bose-Einstein correlations 
may increase (N*^) by ~ 3-10% (see 0). 

The shifts in (N^^) at the had ron level predicted by the models studied thus far 
are given in table 0, where A{N^^) is defined as the change in mean multiplicity 
relative to the 'no reconnection' scenario of each model. From these, it is clear that 
the multiplicities themselves and also the magnitude and sign of the predicted shifts 
are model dependent. 

In this study, the precision with which such tests may be performed is quantified. 
As a starting point for such tests, it was first verified that in the absence of 
reconnection effects {N^^) — "^iN^^^^) in the models pythia and herwig. The 
statistical uncertainty of this test was 0(0.1%). Next, samples of 10^ herwig and 
PYTHIA W+W" events were generated at = 171 GeV with a full simulation of the 
OPAL detector, and realistic event selections were applied for both W^W^^qqqq and 
W+W~^qq^£ {£ = e, n and r). The efficiency in each case was ~80%, while the 
purity is - 80% for W+W-^qqqq and - 88% for the W+W'^qq^^ channel. 

The resulting (uncorrected) charged multiplicity distributions for the hadronic and 
semi-leptonic channels are shown in Figs. |](a) and|l](b), respectively. The simulated 
data correspond to an integrated luminosity of 10 pb~^ at ^/s = 171 GeV, i.e. that 



f In studying these models, no retuning was performed when reconnection was enabled. 
J This was suggested by B R Webber, as a partial emulation of the model of reference [ 
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Table 1. Mean charged multiplicities, (N^^^), and predicted shifts for various models 



model 






A{N^) (%) 


PYTHIA 


normal 


38.64 






type 1 


38.21 


-l.litO.l 




type II 


38.39 


-0.7±0.1 


HERWIG 


normal 


37.07 






reconnected {P = ^) 


37.25 


+0.5±0.1 




reconnected (P = 1) 


38.38 


+3.5±0.1 


ARIADNE 


normal 


38.14 






reconnected 


37.07 


-2.8±0.1 




Figure 1. Inclusive charged multiplicity distributions with 10 pb ^ of fully 
simulated data, with background indicated hatched, at y's = 171 GcV for (a) 
W+W^— >qqqq, and (b) W+W~— >qqfi7^ events, (c) Variation of (N^^) with v^- 
(d) Luminosity dependence of the statistical uncertainty of {N^^) — 2{N^^'^^) (units 
of multiplicity). 



delivered by lep in 1997. In both distributions, the expected background is shown as 
a hatched histogram. The significant level of Z°/7 qq background is apparent in 
the fully hadronic channel. 
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To extract the mean charged muhipUcity at the hadron level at a fixed centre- 
of-mass energy from such distributions, one can apply a simple correction, based on 
Monte Carlo, to the observed mean value, after subtracting the expected background 
contribution. An alternative is to carry out a matrix-based unfolding procedure using 
the event-by-event correlation between the charged multiplicity at the hadron level 
and that observed in the detector after all analysis cuts have been performed. A 
separate correction for the effects of initial state radiation are necessary in this latter 
case. A third alternative is to integrate the fragmentation function but this is not 
discussed here. 

Based on the the simulated data in Fig. |l|(a) and (b), the expected statistical 
uncertainty on the difference {N*^) ~'2{N^^^'^) for an integrated luminosity of 10 pb^^ 
is 2.2 units, or 5.7% on (N^^). The evolution of the precision of such difference 
measurements with more data is estimated using the following assumptions. Firstly, 
the distributions of A^^*,^ and N^^^'^ are seen to be relatively insensitive to changes in 
centre-of-mass energy once away from the threshold region, as illustrated by the energy 
dependence of (N^^) in Fig. |l](c). Therefore both the mean and the corresponding 
rms are assumed constant at their 184 GeV values. Secondly, above = 184 GeV 
the W"'"W~ production cross-section is predicted to vary by less than 10% in the 
region up to ^/s < 200 GeV, and so a constant cross-section of 16 pb is assumed. 
Thirdly, it is assumed that the selection efficiency at 171 GeV may be maintained 
at higher energies. The expected background cross-section is not important as it is 
subtracted in performing the measurement. Given these assumptions, the dependence 
of the expected statistical error on the difference, (5((A^,t> - '^i^X^'^)), is shown as a 
function of integrated luminosity in Fig. |^(d). 

Typically in such multiplicity determinations, systematic effects become significant 
below a statistical precision of 0.5 units of multiplicity. Uncertainty in the modelling 
of 4-jet like Z°/7 qq background with parton shower Monte Carlos in the fully 
hadronic channel may become a significant systematic. 



3. Event shapes 

Global event shape variables have been considered in earlier studies as potential 
signatures for reconnection j2[ |[ ||. In most studies the predicted effects on such 
observables induced by reconnection has been sufficiently small that detection would 
be marginal, even with an integrated luminosity of 500 pb~^. 

The choice of a 'no reconnection' reference sample with which to compare data 
deserves some thought. In trying to find sensitive observables, using the models 
alone is ideal. However, once possible signatures have been developed, and one 
starts to search for effects in data, it will be invaluable to have a well defined 'no 
reconnection' reference sample in data to reduce model and tuning dependence, lep i 
data provide a high statistics reference, but additional assumptions are necessary in 
either extrapolating energy scales, or in combining pairs of —> qq to emulate 
W+W~— >qqqq events without reconnection. It is also necessary to assume that data 
recorded and processed by the detectors before 1996 can be directly compared with 
those recorded near the end of the lep 2 programme. For some signatures, the ideal 
reference data are W"'"W~— >qq£F£ events. However, this sample has only limited size 
and the comparison may require the association of pairs of jets with Ws in the fully 
hadronic channel, a procedure which experimentally introduces more uncertainty. In 
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Figure 2. (a) Effect of typical experimental selection on thrust distribution, and 
(b) hadron level rapidity distribution in ARIADNE for p < 1 GeV. 



the following, all changes are relative to the 'no reconnection' version of each Monte 
Carlo model and all samples are W"'"W~^qqqq. 

This study compares the differences in the rapidity distribution of charged 
particles, dNch/dy^ relative to the thrust axis of each event, in the central region, 
\y\ < 0.5 and for all y, as suggested in ||, |^. As the effects are expected to be 
more pronounced for softer particles, the distribution is studied for three momentum 
ranges, p < 0.5 GeV, p < 1 GeV and all momenta. It has been suggested [||, ^ that 
reconnection effects may be more pronounced in specific topologies where the quarks 
from different Ws are close to one another, therefore events are also studied for all 
thrust values and for T > 0.76. One aspect not considered in previous studies has 
been the effect of applying a realistic event selection, which is necessary to reduce 
the large background {a{Z° — > qq) ~ 20ct(W"'"W~— >qqqq)). As this is dominated 
by two-jet like events, the efficiency for selecting W"*"W^— ^qqqq events in a similar 
configuration is relatively small, as illustrated in Fig. §(a); ~ 38% of W+W^^qqqq 
events selected satisfy T > 0.76, falling to < 0.05% for T > 0.92. 

In the rapidity was studied relative to the axis bisecting the two di-jet axes, 
as a function of the angle separating these axes. Experimentally, without any reliable 
charge identification algorithm to separate quarks from anti-quarks, the specific angle 
proposed in |^ must at best be folded in experimental analyses, and also requires 
pairing of jets into Ws. While the reliability of associating the 'correct' jets together is 
possible with moderate efficiency using kinematic fits, selecting high thrust events was 
used in the current studies for expediency and simplicity. As the shifts in Mw expected 
are modest compared to the experimental mass resolution on an event-by-event basis, 
it is worth considering the use of kinematic fits in which our current knowledge of 
Mw is applied as a constraint, in a similar way to that used by experimental TGC 
analyses. 
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Hadronic events were generated using the models pythia, herwig and ariadne, 
with and without a simulation of the opal detector, and dNch/dy studied within the 
ranges of y, p and T described above. A smearing simulation of the opal detector, 
which is reliable for studies in the W^W^— s-qqqq channel and necessary to achieve 
the relatively high statistics required, was used herein and also to estimate shifts in 
Mw- 

As an example of how the differences may be concentrated in restricted rapidity 
intervals, Fig. ^(b) shows the dNch/dy distribution for p < 1 GeV in ariadne, for 
events with and without reconnection. Changes in charged multiplicity, A(A'^j^), 
within given p and y intervals are summarised in Fig. ^(a) for each of the models 
introduced in table ||, without detector simulation. The left (right) hand side of the 
figure shows the percentage change in (N^^) for the three momentum ranges considered 
for all y {\y\ < 0.5). The leftmost points in this figure correspond to the results of 
table I. Fig. |(b) gives the analogous results for T > 0.76. For illustration, statistical 
errors corresponding to an integrated luminosity of 500 pb~^ are given for the 'herwig 
colour octet' model. 

It is seen that in all models the magnitude of the change increases when only low 
momentum particles are considered. Applying a thrust cut such as T > 0.76 rejects 
^ 40% of events and may change (N^^) by up to two units, but differences relative 
to the 'no reconnection' scenarios are essentially unchanged, therefore the sensitivity 
is reduced. The predicted maximum statistical significance of A{N^^), as well as its 
sign, depends strongly on the model, varying from ^ 6a for ariadne and the herwig 
'colour octet' model, ~ 3.5a for pythia type I, ~ 2a for pythia type II, down to 
~ 0.8cr for HERWIG. The point of maximal sensitivity is indicated (square markers) 
for each model in the figure. Similar trends were observed in studies with detector 
simulation but typically A{N^^) was found to be ~ 50% smaller. 

It may be possible to increase the sensitivity to reconnection effects using charged 
multiplicity based methods, by considering particle distributions relative to the 
W+W~ decay axis, as reconstructed using kinematic fits. In an alternative 
multiplicity signature ('interjet multiplicity') was introduced, having similar sensitivity 
to integrating dNch/dy for |y| < 0.5. This interjet multiplicity was similar in idea 
to methods normally used to quantify the 'string effect' in 3-jet e+e" events. It 
was suggested that this be studied further, using the shape of the particle density 
distribution as a function of the angular separation between jet pairs, rather than 
restricting the study to the integrated particle density in the fixed angular regions. 
However, the 4-jet case is somewhat more complex than the familiar 3-jet case, being 
non-planar, and so this was not pursued during the workshop. 

4. Shifts in Mw 

Extracting Mw from the decay products observed by experiments is non-trivial, 
requiring much attention to bias induced from effects such as initial state radiation, 
detector calibration, imperfect modelling of the underlying physics processes and of 
the apparatus. In comparison to this, estimating a shift which could result from the 
effects of reconnection phenomena is straightforward, as the value of interest is the 
relative shift between Mw determined in two different scenarios of the same model. 
The absolute value of "Mw" obtained is not central to these studies. However, there 
are still many uncertainties inherent in such studies, such as sensitivity of the method 
used to extract Mw to changes in ^/s, to tuning of the Monte Carlo models (e.g. 
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Figure 3. Fractional change in charged multiplicity as a function of maximum 
particle momentum, in two rapidity regions, for (a) all T, and (b) T > 0.76. See text 
for details. 



virtuality cut-offs in the parton shower development), to treatment of combinatorial 
background and ambiguous jet-jet combinations, and the range over which fitting is 
performed to name but a few. 

In these studies, the method used to extract Mw followed closely that used 
by OPAL for its preliminary Myv results using 172 GcV data. In this, events with 
detector simulation are first selected using the same procedure as noted earlier. Four 
jets are formed using the k± jet finder, corrected for double counting of energy 
within the apparatus, and a paramctrisation of the errors on the measured jet 4- 
momenta is carried out. A five-constraint kinematic fit, in which the jet-jet pair 
masses are constrained to be equal, is performed for each of the three possible jet-jet 
pairings, event by event. A mass distribution is constructed using the mass from the 
combination having the highest probability from the kinematic fit in each event if this 
has probability greater than 1%. A second entry is also admitted if the second most 
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Figure 4. Typical mass distribution with fit results with detector simulation and 
event selection. 



probable fit result has probability greater than 1% and within a factor of three of 
the highest probability combination. The aim of this is to include additional mass 
information for events in which the most probable fit combination is incorrect. In 
such events, these two masses are essentially uncorrelated. A typical mass distribution 
formed by this procedure is given in Fig. ^. 

This method was applied to simulated events from each of the models in turn, 
and the shifts obtained are summarised in table where uncertainties on these shifts 
are statistical. The ariadne model predicts a modest shift in mass of approximately 
50 MeV. No significant shift is predicted by the models pythia and herwig. In an 
earlier study, performed in a similar way, significant shifts were determined The 
PYTHIA and ARIADNE models considered in the present study were also included in ^ , 
albeit with different model dependent parameters and looser event selection criteria. 



Table 2. Table of shifts in My/ for each model. 

(AMw> (MeV) 



model 






selected events (e ~ 8 


0%) all events 


PYTHIA 


type I 




+ 18±11 


+11±11 




type II 




-13±11 


-19±11 


HERWIG 


reconnected (P = 


1) 


-16±16 


-19±16 




reconnected (P = 


1) 


+13±15 


+8±14 


ARIADNE 


reconnected 




+51±16 


+51±15 
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One quite plausible explanation proposed was that the difference was due to the 
significantly more stringent event selection currently used. It has been shown that the 
current selection preferentially rejects events having two-jet like characteristics, which 
is where reconnection effects may be expected to be prevalant. The rejection of these 
events does not appear to be the reason for small mass shifts, as a similarly small 
effect is observed when all events are selected, as seen in table |^. 

Many possible sources for the difference were investigated in the context of the 
PYTHiA models. Neither changes in the tuning of pythia/jetset by OPAL^to improve 
the description of lep i data, nor the different centre-of-mass energy {y/s = 175 GeV 
in Q) were found to be significant. The current analysis procedure is slightly different 
to that in j^. However, significant shifts are still found when the current procedure 
is applied to the same simulated events used in Conversely, applying the former 
procedure of |^ to the samples herein does not induce a significantly larger mass shift. 

One apparently significant effect was found to be the choice of mass assigned to 
jets in performing kinematic fits. As discussed in this choice is not unique. In 
the analysis of ||^, the hadronic jets were assumed massless whereas in the current 
studies, the measured jet mass was used. Re-analysing the same simulated events 
of Q but assigning measured masses to the jets reduces the mass shifts estimated, 
e.g. shifts quoted in |] of 130 ± 40 MeV (type I) and 50 ± 40 MeV (type II) become 
70 ± 40 MeV and 30 ± 40 MeV, respectively. For comparison, a sample of 200 000 
fully hadronic type I events were generated at y/s = 175 GeV using identical model 
parameters and program versions, and analysed using the procedure of also using 
measured jet masses. This yielded an estimated shift of 46 ± 16 MeV. It should be 
noted that fluctations due to finite Monte Carlo statistics have to be considered when 
comparing with the results of in which samples sizes for the analoguous studies 
were 50 000 events. 

Comparing the results for mass shifts in table ^ with multiplicity shifts in table |l|, 
it can be seen that any relationship between them is model dependent. Furthermore, 
relatively large shifts in the charged multiplicity do not necessarily lead to a significant 
shift in Afw 

5. Future 

The future for experimental studies of colour reconnection is quite open. There is clear 
model dependence in signatures and mass shifts may be smaller than earlier proposed 
1^ , although there are other models available j^, ^ which were not tested in this study 
from which different conclusions may be drawn. A necessary condition for a model 
to be taken seriously is that it should describe the data, therefore tuning of models 
has to be addressed. With the current statistical precision of lep 2 data, none of the 
models has been put to a stringent test. The effect of background cannot be ignored 
in the W+W^^qqqq channel as it proves difficult to remove. More sophisticated 
selections may be developed, but typically these make use of non-trivial correlations 
between observables, which may be poorly described by the models. A particular 
concern is the description of parton shower Monte Carlos to describe the hard, 4-jet 
like background which is selected. The remaining point of note is that given the model 
dependence inherent to such studies, it is most important to develop signatures which 

t Among these, the cut-off parameter, Qo, was increased from 1.0 GeV in the similar investigation 
of ||, to 1.9 GeV. 
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can be tested taking the 'no reconnection' scenario from data themselves. 
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